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Abstract 

We reexamine the usefulness of fitting blended lightcurve models to microlens- 
ing photometric data. We find agreement with previous workers (e.g. Wozniak & 
^ . Paczynski ) that this is a difficult proposition because of the degeneracy of blend frac- 

Q I tion with other fit parameters. We show that follow-up observations at specific point 

^ I along the lightcurve (peak region and wings) of high magnification events are the most 

helpful in removing degeneracies. We also show that very small errors in the baseline 
magnitude can result in problems in measuring the blend fraction, and study the im- 
portance of non-Gaussian errors in the fit results. The biases and skewness in the 
distribution of the recovered blend fraction is discussed. We also find a new approx- 
imation formula relating the blend fraction and the unblended fit parameters to the 
underlying event duration needed to estimate microlensing optical depth. 
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p i' 1. Introduction 

o ■ 

^ I Gravitational microlensing has become a useful tool in measuring the amount of matter along 

^ ' the line-of-sight to distant stars. Since gravitational lensing depends only on mass, microlensing is 

sensitive to all compact forms of matter independent of their luminosity. Thus measurements out of 
the plane of the Galaxy towards the LMC, SMC, and M31 have given important limits/detections 
of dark matter in the halo ( Aubourg, et al. 1993; Lasserre et al. 2000; Alcock, et al. 1993; 1997a; 
2000; Paulin-Henriksson et al. 2003; de Jong et al. 2004), and measurements towards the Galactic 
bulge give important constraints on the mass and distribution of Galactic stars, including those too 
famt to observe directly (Griest et al. 1991; Paczynski 1991; Han & Gould 2003; Udalski, et al. 
1994; Alcock, et al. 1997b; Afonso, et al. 2003; Popowski, et al. 2005; Sumi, et al. 2005). 

The signal of microlensing is a specific transient magnification of a background source star as 
the lens object passes in front of it, and thus microlensing experiments repeatedly monitor many 
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ordinary stars to find microlensing lightcurves. The probability of microlensing occurring to a 
given star is called the optical depth, r and is of order 10~^ or less for many Milky Way lines-of- 
sight. The smallness of r means that microlensing experiments concentrate on very crowded star 
fields where many hundreds of thousands of stars can be simultaneously imaged. This allows many 
lightcurves to be created simultaneously but also results in blending of the source stars together. 
This blending causes two problems in using the detected microlensing events to infer the optical 
depth. First, since each "source object" may contain the light from many stars, the number of stars 
being monitored is not just the number of objects being photometered. Second, the magnification 
profile of a microlensing event is changed when unlensed light is blended with the lensed light of 
the source star. 

In this paper we revisit the problem of blending in microlensing lightcurves. There are sev- 
eral methods of dealing with the blending problem. Among these are 1) obtaining high resolution 
images from space, which will usually allow separationof the source object into its different com- 
ponents, giving a direct measurement of the fraction of light from the lensed source (Alcock, et al. 
2001), 2) if the unlensed light is not exactly centered on the lensed source, then the centroid of the 
light will shift during the microlensing event allowing limits on blending to be placed (Alard, Mao, 
& Guibert 1995), 3) if the lensed source is a different color than the unlensed light then a color shift 
will occur as the event proceeds, allowing limits on lensed-light fraction to be made (Alard, Mao, 
& Guibert 1995), and 4) for image subtraction lightcurves, the source can in principle be removed 
and this can help break the degeneracy in some cases (Gould & An 2002). 

However, we will not discuss the above methods in this paper but will focus on the fitting 
and interpretation of the photometric data alone; that is, we include the lensed-light fraction as 
a parameter in the microlensing fit and hope to use the shape of the Ughtcurve to recover this 
information. In principle this allows recovery of the actual event duration, and a measurement 
of the amount of blending in the sample of events, allowing corrections to be made in estimating 
the optical depth. A related and popular method is to calculate lensing optical depth using only 
a subsample of very bright source stars (e.g. clump giants). The idea is that very bright source 
stars are less likely to be blended, and when they are blended, should be blended only by a small 
amount. In this case one would like to use the blend fits only to determine whether or not a given 
event is blended. 

Unfortunately, as pointed out previously (e.g. Han 1999; DiStefano & Esin 1995; Wozniak & 
Paczyriski 1997; Alard 1997, etc.) blended fits tend to be quite degenerate. A lightcurve with a 
small lensed-light fraction looks very much like an unblended lightcurve with a smaller maximum 
magnification and a smaller event duration. As pointed out previously, this means that this fitting 
method will be of limited use in many cases. Our study adds strength to the conclusions of previous 
workers, points out several new problems with blend fits, and makes recommendations on how 
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best to proceed with blend fits for those who choose to do them. We will discuss what happens 
when the microlensing event contains signal from other physical effects such as weak parallax or 
binary effects. These effects are not rare, and since the difference between blended and unblended 
lightcurves is small, even an almost undetectable real deviation from the standard point-source- 
point-lens Ughtcurve can render blend fit results meaningless. 

The plan of the paper is as follows: In § 2 we define our notation and discuss the similari- 
ties and differences between unblended microlensing and blended microlensing We also give an 
analytic approximation that gives the underlying event duration and peak magnification from the 
lensed-light fraction and the easily measured apparent event duration and maximum magnification. 

In § 3 we discuss the usefulness of blend fits and compare with earlier work. In § 4 we discuss 
the optimal times to take follow-up data in order to improve recovery of parameters from the blend 
fit. In § 5 we discuss the problem of the baseUne magnitude, and In § 6 we discuss the problem of 
non-Gaussian data and whether the errors returned by fitting programs are reUable. 



2. Degeneracies in blended lightcurves: analytic approximations for event duration and 

A. 



When an isolated lens object crosses close to the line-of-sight of an isolated background source 
star, the source is magnified and a microlensing lightcurve is generated with magnification 

where u is the projected distance between the lens and source in units of the Einstein ring radius, 
is the time to cross the Einstein radius, and t is time, with maximum magnification, ^max, 
occurring at imax- 

The most important parameter is the event duration Ie since the optical depth depends upon 
the sum of efficiency weighted event durations: 

events ^ ' 

where the exposure E is the product of the length in days of the observing program and the number 
of observed stars, and e is the efficiency of detecting an event of duration tE- 



When other sources of light are contained in the same seeing element as the lensed source 
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star, the microlensing lightcurve is altered since only a fraction of the Ught is actually lensed: 

A'{u) = hA{u) - h + 1, (3) 

where /u is the fraction of light that is lensed (a.k.a. the blend fraction, i.e. coming from the source 
star) before the lensing event begins.^ Compared with unblended events, blended photometric 

microlensing lightcurves suffer from a smaller maximum magnification, ^^ax shorter event 
duration t'^, as well as from potential color shifts if the blended light has a different spectrum. 

In fact, a blended lightcurve looks remarkably like an unblended lightcurve with different 
values of tE and ^^ax (e.g. Han 1999; DiStefano & Esin 1995; Wozniak & Paczyriski 1997; Alard 
1997). However, as illustrated in Figure 1, this similarity is not perfect and there are differences 
in the shapes of blended and unblended lightcurves. It is these differences which give rise to 
the hope that information about blending can be extracted by fitting hghtcurves with blending 
parameters. If this similarity were perfect, then there would be no use in fitting blended lightcurves 
to photometric data. Figure 1 shows that the shape differences are typically small, meaning that 
extracting blending information will be difficult. Wozniak and Paczynski 1997 (WP) studied this 
in detail and gave regions of the /u, ^max plane where blended fits were useful and where they were 
not. We return to this subject in § 3, but the qualitative results of WP can be seen from Figure 1, 
where low magnification events show a maximum difference between the blended Ughtcurve and 
the best fit unblended hghtcurve of only 1% or so, while the higher magnification events show more 
substantial differences. 

Previous workers have also given analytic formulas relating the measured (apparent) maxi- 
mum magnification A'^^^, and the apparent event duration t'^ to /n, ^max> and Ie- For example, 
Wozniak and Paczyriski (WP) studied the degeneracies by performing expansions of the above 
equations in the limits of small Umin and large Umm and in these Umits give the formulas relating 
the actual values of and ^max to /u and the measured ^^ax and t'^. For small Umm (large ^max) 
they found u'^^^ «min//ii, and t'^ fntE, while in the hmit of large Umm (^max ~ 1) they 
found !^ Umin/fi/'^, and t'^ ^ fi/'^tE- 

DiStefano & Esin (1995), and Han (1999), and Alard (1997) took a different approach, solving 
equation 3 for A^ax and giving the actual tE in terms of t'^ by requiring that the two different 



' Several terms have been used in different ways in the literature for blend fraction, most commonly ft which cither 
means the fraction of light coming from the lens or the fraction of the light coming from non-lens sources. We introduce 
the new symbol fn to avoid the extant confusion of nomenclature. 
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Fig. 1 . — Four example blended lightcurves (solid) compared with the best fit unblended lightcurve 
(dashed), as well as the difference, A^, between them (blend fit minus unblended fit). The bottom 
labeled time axis is in units of the apparent Einstein Ring crossing time, t'^, that is easily available 
from the data and an unblended fit. However, the extent of the time axis is ±4t^, (labeled on the 
top) where t^; is the underlying event duration used in optical depth estimates. Thus the extent of all 
the time axes is roughly 160 days for a typical microlensing event of duration 20 days. Part (a) has 
values: /u = 0.5, Urain = 0.4, u^j^ = 0.634, and t'j^/tE = 0.754. Part (b) has values: /u = 0.2, 
""min = 0.2, u^ij^ = 0.655, and t'^/ts = 0.475. Part (c) has values: /u = 0.5, Umin = 0.03, 
""min = 0-062, and t'^/tE = 0.594. Part (d) has values: /u = 0.05, -Umin = 0.03, u'^^^ = 0.46, 
andt'^/ifi = 0.144. 
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parameterizations give the same amount of time with A > 1.3416. They found: 




1/2 



^max(HDE) — i^'r, 



1 + /ll)/ flh ^E(HDE) = t'E 



(4) 



where u^j^ = u{A'^^), and ui = u{A{A' = 1.3416)), can be found from the inverse of equa- 
tion 1: 



Noting in Figure 1 that the differences between the blended and unblended lightcurves tend to 
be large in the peak, and that the values of t'^ and A^ax found by fitting, we worried that the 
HDE formula, which assumes equality in the peak, might not be accurate. We also wondered about 
the range of applicability of the WP formulas and so decided to test these formulas. We did this by 
fitting artificial blended lightcurves with an unblended source model and finding the best fit values 
of and ^max- We also fit these Ughtcurves with blended source models and correctly extracted 
the input blend parameters. 

As shown in Figure 2, we found that the WP formulas are not very useful over most of the 
parameter range, and that the HDE equations work well only over a restricted range of param- 
eters. For the WP formulas this is not surprising since they were created only to show that the 
degeneracies exist in certain limits. For relatively large lensed-Ught fraction and for relatively low 
values of ^max the HDE equations give a good estimation of the best fit ^^ax ^'e' small 
lensed-hght fraction or high ^max the estimates of these equation can be far off. 

As expected, it is just where the blended lightcurve shape differs the most from an unblended 
fit that the HDE approximations do not work well. The reason can be seen in Figure 1 , where for 
high magnification events and low lensed-light fraction the blended lightcurve differs strongly in 
the peak area, but not so much in the lightcurve middle rising and falling regions. Thus, the best 
fit unblended lightcurve will allow the actual peak magnification to overshoot and compensate for 
these points by undershooting in the middle regions. Since the HDE formula forces the lightcurves 
to match at the peak and when A'^^ = 1.34, it will overestimate the best fit peak magnitude and 
underestimate the event duration. 

By studying many such examples, one can come up with a formula that does a better job of 
relating the best fit ^^ax ^'e ^'^ ^max» tE, and /u in the parameter ranges where the HDE 
formula does not work well. The points in Figure 2 show the best fit values of ^'„ax ^^'^ ^'e 
/u found by fitting artificial blended lightcurves. The dashed lines show the HDE estimates and 
long dash Unes the WP estimates for t'^/tE- At small values of Amax (< 3) the HDE formulas do 
work very well (better than the new formula) and they should be used. However for ^max > 3 the 
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(A) = (2/Vl - l/A^ - 2)^/2 



(5) 
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HDE formulas do not give accurate estimates. To find a better approximation, one can fit a straight 
line to the data for a given Umin and get a formula which fits well except for very low lensed-light 
fraction. Repeating this procedure for different values of Umin, one discovers that the slopes and 
zero points of the linear fits are also quite linear in u^ain- Thus a simple fitting linear formula that 
covers much of the parameter space can be found. However, if one fits a quadratic for the low /u 
events one can get an even better formula which works very well for /u < 0.3. Thus we find an 
approximation: 



t'E/tE 



^HDE, 

^max-0-9785+0.4150/ii 
O.8I53/11+0. 00021 

^;„3^-0.3618+0.2106/ii 
1.0282/11-0.04433 



t', 



£(HDE) 



/tE, 



ifAr 



< 3; 



if Amax > 3 and A',^ < 10; 



if^r 



> 3 and A'^^ > 10, 



(-1.0946txi„in + 0.9418)/ii + l.Mln^in + 0.0564, ifA^^ > 3 and /u > 0.^) 



FCU, 



if^max > 3 and /u < 0.3, 



where 



(1,/11,/n), u 




(7) 



/ 0.02548 1.0626 

1.1914 7.284 

Y-0.8824 -26.58 

and iimin is found from ^max and equation 5. 



-1.6504\ 
-11.50 
44.28 J 



(8) 



In using this formula, one typically starts with measured values of f/^, A^^^^, and an initial 
guess of ^max and uses different (unknown) values of /u, to find the corresponding underlying 
^max and tE- If the value of ^max found using the new fitting formula is smaller than 3, then one 
should use the HDE formula instead. 



The new fitting formula is shown as the solid line in Figure 2 and does better than HDE or 
WP for Amax > 3. Over the range 0.01 < /u < 1.1, and 3 < A^ax < 70 the new fit formula 
gives a typical error in tE (compared with actually fitting the microlensing lightcurve with a blend 
fit model) of around 3% and a maximum error of 9%. For A^ax the typical error is 4% and the 
maximum error is 12%. The HDE formula can be off by more than 50% in tE and 24% in A^ax in 
this region of parameter space. 
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- Comparison of approximation formulas for relating the underlying microlensing event 
tE and maximum magnification ^max> to the blend fraction, /u, and easily measured 
event duration and maximum magnification, t'^, and A'^^. The circles give the actual 
and unblended results from our Ughtcurve fitting program, while the solid fines show 
approximation formula. The short dashed line shows the HDE approximation, while 
dashed lines show the WP approximations in their two limits. Part (a) is for an actual 



2 {A 

max 

2.65). 



5.07), part (b) shows 



.03 (A 



33), and part (c) shows 



A 



-9- 



In summary, we tested the HDE formula, Wozniak and Paczynski (WP) formulas and Equa- 
tion 6 over a wide range of parameters and found the new fitting function works better than HDE 

for all values of f\\ when ^max > 3 and A'^^^ > 1.34, while the old HDE formula works better 
for low values of ylmax and A'^^^. The WP large v4max formula gives tE within 10% only for large 
/ii(> 0.5), and large A^aax, while the other WP formula is not useful except for ^max -C 1.34. 

Since in microlensing experiments the event durations are found by photometric fitting and 
since the optical depth is proportional to the sum of the fit t£;'s, when making corrections for 
blending it is important to properly relate the lens-light fraction of each event to the underlying 
event duration. 



3. Usefulness of blend fits 

Wozniak and Paczynski (1997) (WP) studied the degeneracy of blend fits and concluded that 
in many cases blended and unblended Ughtcurves cannot be distinguished by photometric fitting. 
They described areas of parameter space where blend fits would be useful and areas where they 
would not. While we think that WP did an accurate and very useful calculation, and we agree with 
their conclusion that blend fits are usually not very useful, we wanted to repeat their analysis for 
several reasons. First, WP did not include the baseline magnitude in their fits, reasoning that since 
many measurements are taken before and after the event, the error in baseline magnitude was not 
significant. In fact, we find that error in the baseline magnitude is one of the most severe problems 
in blend fits. We find that errors even at the few percent level can drastically alter the parameter 
values extracted from the fit. Second, WP considered only evenly spaced observations and we 
wanted to consider whether different follow-up strategies could improve the ability to extract the 
parameters. 

In our studies, we find the error in fit parameters three ways. First we create artificial Ughtcurves 
using the theoretical formula and add Gaussian random noise to each measurement. We perform 
blended and unblended fits on these lightcurves using Minuit (CERN Lib. 2003). Second we 
calculate the error matrix by inverting the Hessian matrix as discussed in Gould (2003). Finally 
to understand the effect of the non-Gaussianity of the errors in real microlensing experiments we 
create artificial lensing lightcurves by adding microlensing signal into actual non-microlensing 
lightcurves obtained by the MACHO collaboration, and then fit these. 

Since the method of calculating the error matrix is closest to what WP did, we first give these 
results. Briefly, we calculate the Hessian matrix (the matrix of second derivatives of the light 
curve residuals with respect to each parameter) then invert it. The square root of the diagonal 
elements of the resulting matrix are then the one sigma errorbars of the parameters. This accounts 
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Fig. 3. — Comparison of our results to the previous results of Wozniak and Paczynski one sigma 
limits on fu for the range of apparent Umin (from non-blend fits). The solid Une is from WP, the 
long dashed line is the 1 — cr/,; Umit for a non-blend fit, and the short dashed line is the value of 

/// for which fu + = 1 for a blend fit. In the region below the long dashed line blending is 
detectable at the one sigma level, above the short dashed line blended events are indistinguishable 
from unblended events, and inbetween the two dashed lines detection is marginal. The region 
where blending is distinguishable can be scaled with a (eqn. 9). 

for correlations in the parameters, but not any nonlinearities. WP used a very similar method, 
but used it to calculate the Ax^ instead of the error bars. In figure 3 we show that our method 
brackets WP's. We show limits calculated as both the one sigma lower Umit on fu for an unblended 
hghtcurve and the value of fu that gives fu + c/„ = 1. 



We note, as WP found, that parameter errors scale linearly with 



(9) 



for N points taken during the peak (defined as lasting AtE)^- Thus our results can be scaled for 
other numbers of observations with different values of a. Thus, we find that our results agree with 
those of WP if we assume the baseUne magnitude is known and take a uniform sampUng. 



^This is true for large enough value of N, for small values of iV < 16 parameter errors increase faster than a. 
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Follow-up Region (t/tg) 

Fig. 4. — The equivalent number of uniform follow-up data points required to improve measure- 
ment of fa as much as a single follow-up observation is plotted as a function of when the single 
follow-up observation is taken. In this case Umin = 0.25, fu = 0.25, and 4 follow-up points are 
added. Times with A?equivaient > 1 are the most effective, while times with A?equivaient < 1 are less 
useful. 

4. Follow-up observations 

Figure 1 shows that the difference between blended and unblended lightcurves is not always 

uniform across the lightcurve. So if one wanted to plan follow-up observations to improve the accu- 
racy of the blend fit, one should concentrate on the regions of the lightcurve where the differences 
are largest. Thus it may be possible to do better than WP suggested with their equal spaced obser- 
vation calculations. To test this hypothesis we calculated the error matrix for blended fits adding 
in follow-up observations at different points on the Ughtcurve. As seen in the Figure 1 examples, 
for any choice of parameters there are five places where the difference Ughtcurves are maximum, 
and therefore where follow-up data is more useful than average: at the peak, in the rising/falling 
portion of the curve, and in the wings. The precise locations change with the choice of parameters 
but for Figure la they are found to be localized near the peak at {\t/tE\ < 0.1), in the falling (or 
rising region) at (0.3 < \t/tE\ < 0.6), and near the baseline at (1.0 < \t/tE\ < 1-5). Observations 
taken between these regions do httle in constraining the parameters. In addition points t greater 
than 2t e are very helpful because they fix the baseUne in our simulated lightcurve. We discuss the 
baseline separately in § 5. In Figure 4 we compare the relative value of added points as a func- 
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Follow-up Region t/tj. 

Fig. 5. — Number of additional uniformly distributed observations required for the same improve- 
ment as 8 focused observations (2 follow-up regions each with 4 observations). The seven con- 
tours are 2 (darkest regions), 4, 6, 8, 10, 12, & 14 (lightest regions). In this case Umin = 0.25 
and fii = .25. Values above (below) 8 indicate an advantage (disadvantage) relative to uniform 
follow-up. 

tion of the time they are added. We find that, in this case, with 40 observations, 4 extra focused 
observations can reduce the error on fu by 7.7%. To get the same reduction of error on fu with 
evenly distributed observations we would need 7 observations, in other words, each added focused 

observation is equivalent to increasing the sampling by 1.75 points over 4 t^. The numerical value 
of the extra effectiveness obtained using focused versus evenly spaced observations varies with 
underlying parameters and the total number of added points. 

Precise follow-up measurements at multiple focused locations can improve the determination 
of even more as they further constrain the shape of the lightcurve. In order to see the effect of 
adding multiple follow-up observations at two distinct times we compare this with adding evenly 
spaced observations. In Figure 5 we plot the increase (or decrease) in effectiveness of extra focused 
observations as a function of the two times at which they are taken. The contours around the Ught 
areas show regions of increased effectiveness, while dark areas show areas where the focused 
observations are less valuable than evenly spaced observations. In this example the effectiveness is 
increased by up to a factor two. 

It is important to note that with more observations or higher accuracy in each follow-up region 
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the advantage per added observation is reduced and the relative values of the various minimums 
vary, though they stay in roughly the same place. For practical use it is important to note that the 
time of the optimum second follow-up observation(s) varies with the time of the first follow-up 
observations(s). In practice one would need to calculate optimum observing times for an event in 
progress as a function of all the previous measurements. 

One problem with the above approach is that without knowing the underlying parameters, 
particularly t^;, it is difficult to predict the best times to take follow-up data. To test if a practical 
experiment could be designed to take advantage of focused follow-up data we simulated an ex- 
periment. First we generated lightcurves with 80 points over 8tE with .05 Gaussian errors at the 
baseline drawing fu randomly from the interval [.01, 1) and Umin randomly from the interval (0, 1) 
requiring ^^ax > 1-34 in the HDE approximation. We also adjusted to keep t'^ ~ lOdays also 
using the HDE approximation. We then generated 9 follow-up observations over 3 days at the peak 
and fit the first half of the light curve plus the follow-up data. From this first fit we calculated 
the optimum times for two more bouts of follow-up. We generated these, both with 9 observa- 
tions over 3 days, and then fit the entire lightcurve with the added 27 points. We also generated 
27 points of follow-up uniformly distributed over the 20 days starting at the peak, added it to the 
initial lightcurve and fit the resulting data. To see the relative improvement for the two methods 
we calculate a parameter ( = (//^^ ^ — /«)/(/// ^ ^ ~ fu)y which is the ratio of the error 
in blend fraction given by focused observations to the error in blend fraction given by uniform 
follow-up sampling. We plot the distribution of ( in Figure 6 finding that our strategy gives an 
improvement < 1) for 71% of the events and a worsening in 29% of the events. We find a sub- 
stantial improvement (C, < .5) for 45% of the events, and and even larger improvement (( < 0.1) 
18% of the time. Thus we conclude that for the the same amount of observing time we can make a 
more accurate measurement of fu by focussing the follow-up observations. 

In summary we find that observations concentrated at a few times can constrain the microlens- 
ing parameters as well as many measurements distributed throughout an event. The best place for 
these measurements are at the peak, in the falling/rising portion, and in the wings with regions 
between where added observations do no good. In most cases it is possible to constrain the event 
parameters well enough with the first half of the data and some follow-up observations near the 
peak to predict the last two optimum observing times. 

5. Baseline Magnitude 

The baseline magnitude of a lightcurve can in principle be very well determined since many 
measurements can be taken before or after the microlensing event. WP assumed that this was the 
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Fig. 6. — The ratio ( showing the advantage of a focused follow-up strategy for 71% of events 
< 1 - unshaded region). For 15% of the events in our simulation are outside the range of this 
plot Id > 2. 



case and so did not include the baseUne magnitude as one of their fit parameters. In real microlens- 

ing surveys, however, it may be that the error in average magnitude is not entirely statistical, and 
may not average down as expected. There may be a systematic limit to the accuracy with which the 
baseline magnitude can be determined. In fact, detectors and telescope systems drift over time and 
so measurements made much later may actually reduce the accuracy of the baseUne magnitude. To 
investigate the importance of the baseUne magnitude, we created artificial Ughtcurves without any 
errors and fit them with a model with a fixed value of baseline magnitude that differed from the 
actual baseline magnitude by various amounts. Our results are shown in Figure 7. We find that 
the dependence on baseline is very strong for low amplification events and not as strong for higher 
amplifications events, but in any case even a 2% error in baseline magnitude determination can 
strongly bias the recovered blend fraction. 

Next, to see how well baseline magnitudes converge in real data, we used the MACHO col- 
laboration database of random stars (Alcock, et al. 2000). We looked at the x^/f^dof of a fit to a 
constant lightcurve for our real Ughtcurves and compared this to simulated ideal Ughtcurves with 
the same number of points and Gaussian errors. For the simulated Gaussian Ughtcurves we find 
the x^/iT'dof distribution peaked near unity and distributed as expected, but for the real data the 
distribution of x^/"do/ is much broader. These two distributions are shown in Figure 8. 
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Fig. 7. — The recovered fu for unblended light curves as a function of an input baseline magnitude 
(fixed at a given value). Forty points over 4 are used. 

As an estimate of the error in the baseline which arises due to the systematic drift and non- 
Gaussian nature of the magnitude errors, we calculated mean and median for the points in each of 
146 lightcurves in MACHO field 119, one of the most frequently observed fields. We found that the 
distribution of mean minus median had a dispersion of 1.3% indicating that the error in the baseline 
flux is ~ 1.3%. Referring to Figure 7 we see that for a typical event with Umin = 0.5, this implies 
a typical spread in fu of 0.18 due to baseUne alone. Since half of all events have Umin < 0-5, half 
of all events will have an even larger bias. For more sparsely sampled fields this dispersion due to 
error in baseUne fit would be even larger. 

6. Errors in Fit Parameters 

From Macho Project data (Table 6 of Popowski et al. 2005) it seems blend fits return biased 
parameters. For the set of Macho clump giant events, which are believed to be minimally blended 
from their positions on the color magnitude diagram, many are best fit with blending. If the events 
are not blended then a systematic bias in the fits must make them appear to be blended. A system- 
atic bias in recovered lensed-light fraction would lead to a bias in the optical depth as well. The 
MACHO collaboration investigated the blending of their clump giant sample and decided to use 
the parameters from the unblended fits. They also used a subsample of events that were less likely 
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Fig. 8. — Actual distribution of 'X^ /ndof (top) for MACHO data and theoretical distribution of 
X^/ndof (bottom) 



to be blended to check for a bias due to blending and found no such bias. 

To test for a systematic bias we generate 1000 lightcurves with Gaussian errors for each of 
three different values for the error on each point: a = 0.01, a = 0.05, and a = 0.15. The 
recovered lensed-light fractions for these events are shown in Figure 9. As the error on individual 
datum increases the distribution of becomes increasingly skewed. We find that while the mean 
fli may not decrease, the most probable value does decrease. This reduction in the mode is at least 
partially compensated by the large tail of the distribution with > 1, but for the small number of 
events a microlensing experiment observes it is unlikely that many of the few events with » 1 
will be observed. Even if one event with //^ » 1 is observed it may be ignored as it is an unphysical 
value of the parameter, thus leading to an underestimate of the average value of fu. Thus we find 
that as the errors in measurement increase blend fitting becomes more and more likely to return 
biased results. The direction of the bias is more often toward small values of Thus events that 
are in reality unblended become more and more likely to return fit values implying that they are 
heavily blended. 
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Fig. 9. — Recovered fu for data with Gaussian errors of 0.01,0.05, and 0.15. As the errors on 
individual data points increase the distribution becomes increasingly skewed with the mode shifting 
toward for larger errors. Also note that 11% of the a = 0.05 events and 24% of the d = 0.15 
events had recovered //^ > 2 while none of the events with a = 0.01 were fit best with > 2. 

7. Conclusions 



We find agreement with previous workers that blend fits are problematic, but can be useful 
especially for high magnification events. When performing blend fits it is helpful to get extra mea- 
surements near the peak and at other specific points along the lightcurve. We find that if care is 
not taken in the treatment of the lightcurve baseline magnitude the fit results can be severely bi- 
ased and in real data the errors returned on fit parameters should be treated with caution. We find 
that blend fits return a biased, skewed distribution of the underlying parameters tending to indi- 
cate more blending than actually exists. Finally, note that when the microlensing event contains 
signal from other physical effects such as weak parallax or binary effects blend fits can yield unre- 
liable results. These effects are not rare, and since the difference between blended and unblended 
lightcurves is small, even an almost undetectable real deviation from the standard point-source- 
point-lens lightcurve can render blend fit results meaningless. 
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